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Abstract. In high energy heavy ion collisions as well as in astrophysical objects like magnetars extreme 
magnetic field strengths are reached. Thus, there exists a need to calculate divers QED processes to all 
orders in the magnetic field. We calculate the vacuum polarization graph in second order of the electric 
field and all orders of the magnetic field resulting in a generalization of the Euler-Heisenberg Lagrangian. 
We perform the calculation in the effective Lagrangian approach of J. Schwinger as well as using modified 
Feynman rules. We find that both approaches give the same results provided that the different finite 
renormalization terms are taken into account. Our results imply that any quantitative explanation of the 
recently proposed Chiral Magnetic Effect has to take 'Strong QED' effects into account, because these 
corrections are huge. 
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1 Motivation and Introduction 

Long time ago J. Schwinger derived in a seminal paper 
[JJ the effective non-linear Lagrangian for constant elec- 
tric and magnetic fields in all orders. Later, the exact 
fermion propagator in a constant magnetic field was de- 
rived based on this work [2J[3] and used to treat e.g. QED 
processes in magnetars [4,5,6,7,8]. The latter is impor- 
tant, because the observed spectra are strongly affected 
by QED effects and the deduction of e.g. the magnetic 
field strength reached in these objects or the structure of 
the accretion column depends crucially on the quality of 
these calculations. 

In heavy ion collisions even far stronger magnetic fields, 
of the order of \B\ ~ (100 MeV) 2 and above, are gen- 
erated for a short time period. Recently, the STAR ex- 
periment at BNL observed correlations between charged 
hadrons which can best be understood, if one assumes 
that topologically non-trivial QCD effects allow for an ef- 
fective, naively CP-odd coupling of the type E ■ B, see 
[MTU] . This hypothetical mechanism is called Chiral Mag- 
netic Effect (CME). In [llj it was argued that the electric 
field induced by any such effect should minimize the en- 
ergy and thus is determined by the linear term from the 
coupling to G°„G aMI/ and the field energy term, quadratic 
in E. Therefore, whatever the precise nature of the CME 
might be, it will be affected strongly by pure QED effects 
which drastically change the electromagnetic field energy 
for such strong fields. Thus they will alter the magnitude 
of the induced electric field strength and thus the size of 
the charged particle correlations. These QED effects are 
even important for all charge correlations, independently 
of whether the CME is confirmed by future measurements 



or not. 

In principle these calculations should take the detailed dy- 
namics of heavy-ion collisions into account. In the present 
work, however, we only discuss the case of constant fields, 
which already involves some conceptual problems. 
In the CME the induced electric fields are relatively weak, 
roughly of the order of (10 — 20 MeV) 2 , such that it is suf- 
ficient to determine the contribution which is quadratic in 
E but includes all orders in B. 

The main technical problem we are facing is the follow- 
ing: Schwinger's elegant calculation is based on his proper 
time formalism, which from the very beginning expresses 
the effective Lagrangian in terms of the gauge invariant 
field strength tensor F^ v ', For many dynamical applica- 
tions one does need, however, the fermion propagator for 
which an explicit form is given e.g. in [3J. We, therefore, 
reproduced the Schwinger result also in that formalism, 
which was highly non-trivial and actually required a more 
careful definition of the exact fermion propagator than the 
one given in [3J. 

These technical problems are also reflected by the litera- 
ture. In many papers the effective action result of Schwinger 
is used to analyze specific situation beyond the weak field 
limit, e.g. [12]. In others the problem is discussed that its 
expansion in powers of the fields leads to an asymptotic se- 
ries, which is very difficult to handle and discouraged some 
applications. However, it seems that this is not a funda- 
mental problem but just an unlucky choice of expansion 
as was shown e.g. in [T3J. In this paper we will therefore 
avoid any expansion. This will lead us to expressions con- 
taining the ip function, which indeed has an asymptotic 
expansion involving Bernoulli numbers but has perfectly 
reasonable properties if not expanded. 
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In e.g. |14| it was discussed with great clarity that finite 
regularisation terms have to be treated with care to avoid 
misinterpretations, but that they do not pose a problem 
of principle. In our calculation we are interested in higher 
order terms for which such problems do not occur and 
adopt for the leading terms just the standard results. 
The outline of this paper is as follows: In section 2 we 
shall present the effective Lagrangian calculation within 
the Schwinger formalism and in section 3 we discuss the 
problems encountered when trying to do the same calcu- 
lations with Feynman rules. In section 4 we will conclude 
and discuss our findings in the context of the CME. 



2 The calculation using Schwinger's 
effective Lagrangian formalism 

Schwinger's expression for the effective Lagrangian of con- 
stant electromagnetic fields reads pQ 
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= -— r / dss- 3 exp(-m 2 s) 
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where 



X = (2{F + ig))* 
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T=\{B 2 -E 2 ) 



(2) 
(3) 

(4) 



and m is the mass of the considered Dirac field, here of 
the electron. With 
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the real and imaginary parts read 
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and Eq.((T]) simplifies to 
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For renormalisation one has to subtract the logarithmic 
divergence and one has to decide on the finite renormal- 
isation one chooses. In principle one could e.g. subtract 
the contribution for any fixed magnetic field Bq. However, 
there is no good reason to introduce such an additional pa- 
rameter. We follow Schwinger in subtracting the limiting 
case for vanishing magnetic field 



C 



B=0 



1 

8^ jo 

4\ 



dss 3 exp(— m 2 s) 



e 2 s 2 



B 2 - E A 



0(£ 4 ) 



(11) 



and thus obtain for the effective, renormalised Lagrangian 
of second order in E and all orders in B 
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The sign convention was chosen such that the contribution 
of e.g. V e ff to the energy density is 



E 



dE 
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The behaviour of the potential V e ff + Vfff i s shown in 
Fig. [T] and Fig. [5] It is easy to verify that for small B 
fields one reproduces all contributions of second order in 
E to the Euler-Heisenberg Lagrangian 
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Fig. 1. Potential energy correction in leading order of the 
electric field as a function of the angle between the electric and 
magnetic field and fi = m 2 /(eB); for better visualization of 
the behaviour at fi — in this plot the potential is multiplied by 
fi, and additionally normalised to the corresponding prefactor 
of the Euler Heisenberg lagrangian e 2 E 2 / (72tt 2 ) . 
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Fig. 3. Relative magnitudes of relevant terms in the lagrangian 
which are quadratic in the electric field as a function of fi — 
m 2 /(eB); fj, is in heavy ion collisions of the order of to 
1(T 5 . 
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with a numerical factor \n which is given in Table[T]for N 
up to 8. It decreases approximately exponentially with the 
order of the expansion N. To get the corresponding terms 
in the energy density one has to multiply the V e ff,N by 
an additional factor of 2N — 1, of course, due to Ea. lfTS]) . 
The result can be tested again using the Euler-Heisenberg 
Lagrangian. For N = 2 and small magnetic fields it repro- 
duces exactly the missing term of order E A in Eq. (flU|) 



Fig. 2. Potential energy correction in leading order of the elec- 
tric field as a function of fi — m 2 /(eB) for values of the angle 
between the electric and magnetic field = and O = it/2; 
in this plot the potential is normalised to the corresponding 
prefactor of the Euler Heisenberg lagrangian e 2 E 2 /(72 7T ). 
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For the CME we are most interested in the case = 0, 
i.e. in V e //(-B, E). The extension of £reL to higher orders 
in E is quite straight forward. For the term of interest one 
gets to arbitrary order N > 1 
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In heavy ion collisions B is much larger than m e such 
that one can use the asymptotic expansion of the "0 func- 
tion 
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For B = (lOOMeV) 2 one finds V e ff(B,E) ~ 4.5E 2 , nine 
times the free electric field energy. Plots of the relevant 
quantities are shown in Fig. [3] As argued in [TT] this in- 
crease of the total electric field energy by a factor 10 will 
in turn reduce the induced electric field strength substan- 
tially. Even more important is the fact that V e // grows 
linearly in B, just as the driving topological term. If it 
had grown faster (e.g. like B 2 as for the Euler-Heisenberg 
Lagrangian), the CME would have died out with increas- 
ing B. As it is the induced electric field is only little B 
dependent for strong B fields. However, one should keep 
in mind that our result applies only to the situation of 
constant fields, which is not a good approximation for a 
heavy ion collision. 
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Table 1. Coefficients of the expansion to higher orders in the electric field 
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4 


nan 


45 


945 


4725 


93555 


638512875 


18243225 


162820783125 



3 An alternative calculation using the 
explicit form of the fermion propagator 

Schwinger's approach is the most efficient one, when the 
aim is to derive the effective higher-order photon action. 
However, one often wants to study dynamical quantum 
processes in a constant background field. To do so, one 
needs the modified Feynman rules in such a background 
field. In this section we will try to re-derive the result 
just obtained for (9 = with the Schwinger approach, 
namely the effective action in second order in the electric 
field and all orders in the magnetic one, using the explicit 
form of the fermion Feynman propagator given in [3] . This 
exercise is meant to demonstrate that one can actually do 
so, but also that special care is needed with respect to 
finite renormalisation terms. The propagator reads 

eB 



be linked to the existence of Landau orbitals in a constant 
magnetic field. The Fourier transform of Eq. (|23l) gives typ- 
ically Gaussians of the form: 
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where an e prescription was introduced such as to make 
a Wick rotation in /i possible, which lies at the heart of 
this approach. It will always make the results finite and 
thus includes implicitly a regularisation. As mentioned in 
the Introduction different renormalisation schemes differ 
by finite renormalisation terms. In the Geprags et al. ap- 
proach this ambiguity might show up as an ambiguity in 
the choice of the e prescription, but this was not explored 
so far. To simplify the notation the e will be suppressed 
from now on. 

The situation is rather complicated. The e lcot ~ function 
has a countably infinite number of essential singularities 
on the positive real axis at 



z N 
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, ViV e N. 
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Thus the analytic continuation is highly non-trivial and we 
will show that for our case this prescription leads to a dif- 
ferent finite renormalisation than the Schwinger formula. 
(The leading logarithmic divergence is naturally renor- 
malised in the same way.) Physics- wise the problem can 
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and for e.g. tan(e_B/2/x) = tt one gets identical weight 
factors for all Landau orbits with 
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with the anomalous magnetic moment a/. This can also 
be interpreted more intuitively as follows: UV regularisa- 
tion is concerned with the short distance behavior. For 
arbitrarily small eB > all classical Landau orbits return 
to the starting point and thus give an unsuppressed con- 
tribution. In contrast Schwinger subtracts only the B = 
contribution. One cannot expect that the \i integrals 
from Ea. (|23p . the sum over Landau orbits and the limit 
B — » commute. The Geprags et al. result should still 
give the correct answer, but only up to finite renormal- 
ization terms. Therefore, we will proceed as follows: we 
will apply this approach disregarding the problems just 
discussed and will adjust the finite regularisation terms 
in the final expression to the Schwinger result. (We were 
not able to find an e prescription which would automati- 
cally result in an expression in the same renormalisation 
scheme as Schwinger's approach.) Luckily, the changes to 
be made are rather obvious. In any case, we find it very 
reassuring that up to this finite regularisation term both 
approaches lead to the same result. 

We treat the magnetic field B exactly, i.e. as part of 
Ho and the electric field as perturbation, i.e. as part of 
Hint and equate the second order of the expectation value 



(B\T{ 



,ifd\C int ( X ) 
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with the first order of the effective electromagnetic La- 
grangian we are seeking (which is of second order in the 
electric field E), i.e. 



£e//(Xl 



2 



rf 4 X2Ant(Xl)Ant(X2) > \B) 



(30) 

We basically calculate the vacuum polarization graph de- 
picted in Fig. @] for a constant E field in the background 
of an external constant B field. For simplicity we only dis- 
cuss the case E = Ee Zl B = Be z most relevant for the 
CME. 
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^§(X',X) 




S(x,x7 

Fig. 4. Sketch of the vacuum polarization Feynman diagram. 



For constant fields, rather than point photons, the propa- 
gator from xi to X2 fulfills 



S(xi,X2) = S 4 (\i - X2) 



(31) 



with the gauge invariant derivative rather than the usual 
one. This introduces an overall gauge factor, which in turn 
leads to covariant derivatives of the form 



eQ 



■F^v {x X ) 



(32) 



see [15] . In [16] this was actually already discussed in con- 
nection with the CME. 

We chose the 4-coordinates (x\) '■= (t,x,y,z), (X2) '■= 
(t 1 , x' ,y' , z'). The expression to be calculated is 



= 2 / d ^ 2 E ^ z ~ z') 2 n m ( X i,X2) (33) 



n y.v 



-e 2 Tr(^S(xi,X2)rS( X 2,Xi)), 



(34) 



where we adopted the gauge (A^)(x) ■= (zE, 0,0, 0). 
With the abbreviations 
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Note that the first exponentials of the propagators can- 
celled in Eq. 43E]) so that 77 00 (x', x) = n 0Q (x' - x)- 
Evaluating the exponentials in the trace yields 



cos v — ia 12 sin v. 
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The evaluation of the trace is straight forward and gives 
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with the notation Xi — X2 = (-^j Ax, Ay, Az) = A and 
X = (xi + X2)/2. Now we rotate the /1 and /i integrals to 
the positive imaginary axis, and perform a Wick-rotation 
for At, yielding 



eE\ 2 f°° , f 00 v v 



£eff = TZa / d V / 



47r 2 J Jo *~ Jo sinh sinh £ 

x e 2 V f A / 

X y rf4 ziziz 2 e -^(coth 1 .+coth 1 >)(zi a; 2 +/ly 2 ) 

x e - <<£+<£> (^ 2 + 4t 2 ) 

x (cosh(i> + v)pjl(Az 2 - At 2 ) 



e 4-Y (Ax 2 + Ay 2 ) - -. - . 
2 / smh sinh z/ 



(41) 



After performing all Gaussian integrals this simplifies to 
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which is an even function in B, allowing us to substitute 
H, p, by p — \v\, p = \i>\. In terms of these variables we 
introduce a = p + p and 5 = p — p and perform the 5 
integration to obtain 
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which differs from the result obtained in the last section 
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only in the subtracted renormalisation term. While in the 
Schwinger treatment the latter is B independent 
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the former is not. We correct this by simply adding the 
relative finite renormalisation term 



AC 



Schwinger 
eff 



e 2 E 2 
24tt 2 



due 
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sinh 
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We conclude that the approach from Geprags et al. can 
be used equally well, if the renormalisation for all loop 
graphs is adapted to the usual conventions. 



4 Conclusions 

In this publication we have calculated higher-order terms 
of the effective electromagnetic Lagrangian. While strong 
QED is an interesting and active field in its own right, 
with applications in e.g. the astrophysics of magnetars, 
these studies were specifically motivated by the Chiral 
Magnetic Effects (CME), possibly observed in high-energy 
heavy ion collisions. The postulated mechanism is that 
the extremely strong magnetic fields present in the early 
phase of such a collision in combination with topological 
tunneling in QCD could induce an electric field, subse- 
quently generating specific charge correlations. For any 
such mechanism strong QED effects are so large that they 
have to be taken into account for any quantitative descrip- 
tion. This was demonstrated in our paper for the case of 
constant electric and magnetic fields. The observed charge 
correlations are small, hinting to an electric field strength 
of the order 100 MeV 2 , while the magnetic field strength 
is of the order 10 4 MeV 2 , such that the most relevant term 
is of second order in E and all orders in B. We calculated 
this contribution. 

The result we got implies that strong QED effects increase 
the energy density associated with an electric field by an 
order of magnitude for B — (100 MeV) 2 , thus strongly 
suppressing the size of such fields as compared to naive 
expectations. This result is obviously most relevant for 
the CME, although the approximation of constant fields 
is probably a bad one for the heavy-ion setting where the 
fields change on time and distance scales of several fm, 
which is also the radius of Landau orbits for the B fields 
considered. If the CME is confirmed by future experiments 
at RHIC and especially LHC, where the fields will still be 
stronger and even more Lorentz contracted, one will have 
to develop techniques to treat also the dynamics of the 
strong QED effects reliably. 

For calculations of e.g. Compton scattering in the accre- 
tion column of a magnetar one needs the Feynman rules 
for a magnetic background field, as given e.g. by Geprags 
et al. [3J. Therefore, we performed our calculation also 
with these explicit Feynman rules and showed that the re- 
sult agrees with that calculated in the Schwinger approach 
up to a finite renormalisation term. We see this as an im- 
portant check for the Geprags et al. approach, which also 
illustrates nicely some of the technical problems encoun- 
tered by any such calculation. 

Finally, in the Schwinger approach we generalized the re- 
sult to higher orders in the E field strength. 
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